Multiscaling in Infinite Dimensional Collision Processes 
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We study relaxation properties of two-body collisions in infinite spatial dimension. We sliow that 
this process exhibits multiscaling asymptotic behavior as the underlying distribution is characterized 
by an infinite set of nontrivial exponents. These nonequilibrium relaxation characteristics are found 
to be closely related to the steady state properties of the system. 
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Our understanding of the statistical mechanics of 
nonequilibrium systems remains incomplete, in sharp 
contrast with their equilibrium counterpart. The rich 
phenomenology associated with dynamics of far from 
equilibrium interacting particle systems exposes the lack 
of a unifying theoretical framework. Simple tractable mi- 
croscopic models can therefore help us gain insight and 
better the description of nonequilibrium dynamics. 

In this study, we focus on the nonequilibrium relax- 
ation of an infinite particle system interacting via two 
body collisions. We find that a hierarchy of scales under- 
lies the relaxation. In particular, we devise an extremely 
simple system which exhibits multiscaling in infinite di- 
mension, while in finite dimensions simple scaling behav- 
ior is restored. Furthermore, we show that this behavior 
extends to a broader class of collision processes. 

We are interested in modeling collision processes in 
a structureless, infinite dimensional space. Therefore, 
we place an infinite number of identical particles on the 
nodes of a completely connected graph. Particles are 
characterized by a single parameter, their velocity v. 
Two-body collisions are realized by choosing two par- 
ticles at random and changing their velocities according 
to (ui,U2) (vi,V2) with 
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and < 7 < 1. In other words, the post-collision veloci- 
ties are given by a linear combination of the pre-coUision 
velocities. Both the total momentum {ui + U2 = vi + V2) 
and the total number of particles are conserved by this 
process. In fact, the collision rule (|^) is the most general 
linear combination which obeys momentum conservation 
and Galilean invariance, i.e., invariance under velocity 
translation v v — vq. 

Our motivation for studying this problem is inelastic 
collisions in one-dimensional granular gases 0-|3). While 
the two problems involve different collision rates, they 
share the same trivial final state where all velocities van- 
ish, P{v, t) — > 5{v) when t 00 (without loss of general- 
ity, the average velocity was set to zero by invoking the 
transformation v — > v — (v)). We chose to describe this 
work in sightly more general terms since closely related 
dynamics were used in different contexts including voting 
systems ^,|[ , asset exchange processes |^ , combinatorial 



processes 0, traffic flows [||, and force fluctuations in 
bead packs We will show that multiscaling charac- 
terizes fluctuations in some of these problems as well. 

Velocity fluctuations may be obtained via the probabil- 
ity distribution function P{v, t) which evolves according 
to the following master equation 



dP{v,t) 
dt 
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duidu2P{ui,t)P{u2,t) (2) 



X [S{v - 7U1 - (1 - 7)^2) - S{v - U2)] ■ 

The 5— functions on the right-hand side reflect the colli- 
sion rule and guarantee conservation of the number 
of particles, J dvP{v,t) = 1, and the total momentum 
/ dvvP{v,t) = 0. Eq. (I) can be simplified by eliminat- 
ing one of the integrations 



dP{v,t) 
di 



+ P{v,t) 



^jypiu,t)p[^,t 
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Further simplification may be achieved via the Fourier 
transform P{k, t) = J dv e^^"" P{v, t) which obeys 



dt 



P{k, t) + P(fc, t) = P[7fc, t] P[(l - 7)fc, t]. (4) 



Although the integration is eliminated, this compact 
equation is still challenging as the nonlinear term be- 
comes nonlocal. 

Velocity fluctuations can be quantified us- 
ing the moments of the velocity distribution, 
Mn{t) = J dvv"'P{v,t). The moments obey a closed 
and recursive set of the ordinary differential equations. 
The corresponding equations can be derived by insert- 
ing the expansion -P(M) = E„ ^^Mn{t) into Eq. (|) 
or directly from Eq. (y). The first few moments evolve 
according to Mq = Mi = 0, and 



M2 = 
M3 = 
M4 = 

with the coefficients 
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i„(7) = l-(l-7r-7", 
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and 024 = 67^(1 — 7)^. Integrating these rate equations 
yields Mq = 1, Mi = and 

M2(t) =M2(0)e-''^* 

M3(t) =M3(0)e-'^^* (7) 
The asymptotic behavior of the first few moments sug- 

1 /2 

gests that knowledge of the RMS fluctuation v* = M2 
is not sufficient to characterize higher order moments 
since Mg^'^/w*, M^^'^/v* — > 00, as t — > cx). 

This observation extends to higher order moments as 
well. In general, the moments evolve according to 

Mn + anMn - (^i)'^'"^^ ~ 7)""™Af™M„_™. (8) 

m=2 ^ ^ 

Note that for < 7 < 1, the coefficients a„ 
satisfy a„ < am + an-m when 1 < m < n — 1. 
This inequality can be shown by introducing 
G(7) = a,„(7) + a„_m(7) - an(7) which satisfies G{0) = 
and G{j) — G(l — 7). Therefore, one needs to show 
that G' (j) ^ mlbrn - b„] + {n - m)[bn-m - K] > for 
< 7 < 1/2 with bn = bnil) = (1 - 7)""^ - 7""^ One 
can verify that the &„'s decrease monotonically with in- 
creasing n, bn > bn+i for n > 2, therefore proving the 
desired inequality. Since moments decay exponentially, 
this inequality shows that the right hand side in the above 
equation is negligible asymptotically. Thus, the leading 
asymptotic behavior for all n > is Af„ ~ exp(— a„t). 
Since the a„'s increase monotonically, a„ < a„_|_i, the 
moments decrease monotonically in the long time limit, 
Mn > Mn+1- Furthermore, in terms of the second mo- 
ment one has 

1 _ d _ ^-j" _ 

M„ cc M^- , a„ = ) ^> \ . (9) 
1 - (1 -7) - T 

While the prefactors depend on the details of the initial 
distribution, the scaling exponents are universal. There- 
fore, the velocity distribution does not follow a naive scal- 
ing form P{v,t) — :p^P{^). Such a distribution would 
imply the linear exponents an = at — n/2. Instead, the 
actual behavior is given by Eq. (O) with the exponents 
a„ reflecting a multiscaling asymptotic behavior with a 
nontrivial (non-linear) dependence on the index n. For 
instance, the high order exponents saturate, a„ — > 
for n 00, instead of diverging. One may quantify the 
deviation from ordinary scaling via a properly normalized 

set of indices /3„ — an/a^ defined from M^" {v*)^". 
By evaluating the 7 = 1/2 case where multiscaling is 
most pronounced, a bound can be obtained for these in- 
dices: 7/8, 31/48 < /3„ < 1 for n = 4, 6 respectively. Fur- 
thermore, /3„ ^ 1 — ^"~'^ 7 when 7—^0 indicating that 
the deviation from ordinary scaling vanishes for weakly 
inelastic collisions. Thus, the multiscaling behavior can 
be quite subtle |p^ . 



The above shows that a hierarchy of scales underlies 
fiuctuations in the velocity. In parallel, a hierarchy of 
diverging time scales characterizes velocity fluctuations 

M^" ^ exp(-t/r„), r„ = — . (10) 

These time scales diverge for large n according to 
T„ ~ rt. Large moments reflect the large velocity 
tail of a distribution. Indeed, the distribution of ex- 
tremely large velocities is dominated by persistent parti- 
cles which experienced no collisions up time t. The prob- 
ability for such events decays exponentially with time 
P(u, i) ~ P(u, 0) exp(— i) for u 1 (alternatively, this 
behavior emerges from Eq. (|) since the gain term is 
negligible for the tail and hence P + P = 0). This 
decay is consistent with the large order moment decay 
Mn ~ exp(— t) when n — + 00. 
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FIG. 1. Development of a singularity for a compact initial 
distribution. Shown is the probability distribution obtained 
by simulating the collision process of Eq. (|^) with 7 = 1/2. 
The data represents an average over 200 independent realiza- 
tion in a system with 10^ particles, starting from a uniform 
distribution in the range [—1 : 1]. 

Although the leading asymptotic behavior of the mo- 
ments was established, understanding the entire distri- 
bution P(y, t) remains a challenge. Simulations of the 
7 = 1/2 process reveal an interesting structure for com- 
pact distributions. Starting from a uniform velocity dis- 
tribution, Pa{v) — 1/2 for —1 < v < 1, the distri- 
bution loses analyticity at v — ±1/2. Our analysis of 
Eq. (HI) shows that such a singularity should indeed de- 
velop at w = ±1/2 and it additionally implies the appear- 
ance of (progressively weaker and weaker) singularities at 
V = ±1/4, etc. More generally, for an arbitrary compact 
initial distribution and an arbitrary 7, the distribution 
P{v, t) loses analyticity for t > and develops an infinite 
(countable) set of singularities whose locations depend 
on the arithmetic nature of 7 (e.g., it is very different 
for rational and irrational 7's). On the other hand, un- 
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bounded distributions do not develop such singularities, 
and therefore, the loss of analyticity is not necessarily 
responsible for the multiscaling behavior. 

Asymptotically, our system reaches a trivial steady 
state P{v,t — oo) — S(v). To examine the relation be- 
tween dynamics and statics, a non-trivial steady state 
can be generated by considering the driven version of 
our model External forcing balances dissipa- 

tion due to collisions and therefore results in a non- 
trivial nonequilibrium steady state. Specifically, we as- 
sume that in addition to changes due to collisions, ve- 
locities may also change due to an external forcing: 
^Ihoat = We assume standard uncorrelated white 
noise {^i{t)^j{t')) = 2DSijS{t — t') with a zero average 
{^j) = 0. The left hand side of the master equation (||) 
should therefore be modified by the diffusion term 



dP{v,t) dP{v,t) d^P{v,t) 



dt 



dt 
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(11) 



Of course, the addition of the diffusive term does not 
alter conservation of the total particle number and the 
total momentum, and one can safely work in a reference 
frame moving with the center of mass velocity. 

We restrict our attention to the steady state, obtained 
by setting the time derivative to zero. The corresponding 
Fourier transform Poo{k) = P{k,t = oo) satisfies 



(1 + Dk^)P^{k) = P^[jk] P^iil - 7)fc]. 



(12) 



The solution to this functional equation which obeys the 
conservation laws Poo(O) = 1 and {v) — P^{0) = is 
found recursively 



^'oo(fc) - n n [i + ^'■''(1 - ^)"^'-'^De 

i=0 j=0 
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(13) 



To simplify this double product we take the logarithm 
and transform it as follows 



I + -f^i (l - ^)^<^'-3^ Dk^ 



lnPo.,(fc)^-^5:r. In 
1=0 j=o ^-^^ 

i=0 j=0 ^"'^ n=l ^ 
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n=l i=0 j=0 ^-^^ 

^^^^^^^,„^^^_^^,„^._ (14) 

n=l i=0 

The second identity follows from the series expansion 
ln(l + q) = —J2n>i^^^i~^)"' ^^'^ ^^'^ forth from the 
binomial identity X]j=o Q)-P"''?* = (p + 9)*- Finally, us- 
ing the geometric series (1 — x)^^ — J2n>o -^"j Fourier 
transform at the steady state is found 



-Poo (A;) = exp ■ 



.n=l 



{-Dk'^Y 
na2n{l) 



(15) 



with a„(7) given by Eq. (^. The nth cumulant of the 
steady state distribution k„ can be readily found from 
-Km- Therefore, the odd cumulants 



lnP^{k)=Y.m^ 
vanish while the even cumulants are simply proportional 

to the time scales characterizing the exponential relax- 
ation of the corresponding moments: 



K2n 



{2n-l)\ 



(16) 



Of course, the moments can be constructed from these 
cumulants. Interestingly, a direct correspondence be- 
tween the steady state characteristics and the nonequi- 
librium relaxation time scales is established via the cu- 
mulants of the probability distribution. 

None of the (even) cumulants vanish, thereby reflecting 
significant deviations from a Gaussian distribution. Nev- 
ertheless, for sufficiently large velocities, one may con- 
centrate on the small wave number behavior. Using the 
inverse Fourier transform of (|l^) one finds the tail of the 
distribution 



Poo{v) 



a2 
AttD 



exp 



a2V 



V » VD/a2. (17) 



This in particular implies the large moment behavior 
Man ^ (2n - l)!!(4D/a2)" as n ^ oo. 

To examine how general is the above behavior, we 
briefly discuss a few generalizations and extensions of the 
basic model. Relaxing Galilean invariance, the most gen- 
eral momentum conserving collision rule is 



71 
1 ~ 71 



1-72 

72 



(18) 



Following the same steps that led to (^ shows that when 
7ii72 7^ 0, 1 and when Mi = this process also exhibits 
multiscaling with the exponents a„ = an/a2, where 
an(7i:72) = ^[anili) + On (72)]- When 71 = 1 - 72 = 7 
one recovers the model introduced by Melzak [Q, and 
when 71 = 72 = 7 one recovers inelastic collisions. Since 
o-ni'y) — an(l— 7) both models have identical multiscaling 
exponents. Furthermore, a multiscaling behavior with 
the very same exponents Q:„{'y) is also found for the fol- 
lowing process (ui,M2) (ui—jui^vi+jui) investigated 
in the context of asset distributions Q and headway dis- 
tributions in traffic flows 

One can also consider stochastic rather than deter- 
ministic collision processes by assuming that the colli- 
sion ([1^ ) occurs with probability density cri(7i,72)- Our 
findings extend to this model as well and the multiscal- 
ing exponents are given by the same general expression 
an = a„/a2 with an = J dji J c;72cr(7i,72)a„(7i,72). In 
particular, for completely random inelastic collisions, i.e.. 
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Tt-1 

n+1 



and hence 



(7 = 1 and 7i = 72 = 7, one finds a„ 

n+l 

So far, we discussed only two-body interactions. We 
therefore consider A^-body interactions where a collision 
is symbolized by (ui, . . . , un) (^i, • ■ • , v^). We con- 
sider a generalization of the 7=5 two-body case where 
the post-collision velocities are all equal. Momentum 
conservation implies Vi = u = '^Ui. The master 
equation is a straightforward generalization of the two- 
body case and we merely quote the moment equations 
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M„, ■•■M„„ (19) 



with a„ 



= 1 — A^^ Using the inequality 

o-n-m for all 1 < TO < n — 1, and its kin like 



a 



m2 



a 



71 — mi —7712 



for all 1 < mi, TO2 < n — 1, 
etc., we find that the right-hand side of the above 
equation remains asymptotically negligible. Therefore, 
Mr, ~ e-°"* and 



M„ 



1 - AT- 



(20) 



Thus, this iV-body "averaging" process exhibits multi- 
scaling asymptotic behavior as well. 

Thus far, we considered the behavior on a mean field 
level, i.e., in an infinite dimensional space. It is natural 
to consider the finite-dimensional counterpart. Specifi- 
cally, we assume that particles reside on a d dimensional 
lattice and that only nearest neighbors interact. Here, 
the above dynamics is essentially equivalent to a diffu- 
sion process p^ . As a result, the underlying correlation 
length is diffusive, L{t) ^ t^/^. Within this correlation 
length the velocities are "well mixed" and momentum 
conservation therefore implies that v ~ L"**/^ ~ t~'^^^. 
Indeed, the infinite dimension limit is consistent with the 
above exponential decay. Furthermore, an exact solution 
for moments of arbitrary order is possible We do 

not detail it here and simply quote that ordinary scaling 
is restored M„ ^ i.e. a„ = a* = n/2. Thus, spa- 

tial correlations counter the mechanism responsible for 
multiscaling. 

In summary, we have investigated inelastic collision 
processes in infinite dimension. We have shown that 
such systems are characterized by multiscaling, or equiv- 
alently by an infinite hierarchy of diverging time scales. 
Multiscaling holds for several generalizations of the ba- 
sic model including stochastic collision models and even 
processes which do not obey Galilean invariance. In this 
latter case, however, multiscaling is restricted to situa- 
tions with zero total momentum. This perhaps explains 
why multiscaling asymptotic behavior was overlooked in 
previous studies MM- Another explanation is that this 



behavior may be difficult to detect from numerical sim- 
ulations. Indeed, in other problems such as multidimen- 
sional fragmentation and in fluid turbulence, low or- 
der moments deviate only slightly from the normal scal- 
ing expectation. 

There are a number of extensions of this work which 
are worth pursuing. We have started with a simplified 
model of a ID granular gas with a velocity independent 
collision rate. One possibility is to approximate the colli- 
sion rate with the RMS velocity fluctuation. This leads to 
the algebraic decay M„ with given by Eq. (||) 

and in particular, Haff's cooling law T = M2 ~ is re- 
covered [0. Our numerical studies indicate that when 
velocity dependent collision rates are implemented, or- 
dinary scaling behavior is restored. One may also use 
this model as an approximation for inelastic collisions in 
higher dimensions as well, following the Maxwell approx- 
imation in kinetic theory 
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